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ABSTRACT 

A unified  theory  of  noise  generation  and  amplification  by 
turbulent  combustion  of  premixed  fuel  and  liquid  fuel  droplets 
has  been  developed  within  the  framework  of  the  fluid  mechanics 
of  the  reacting  gas.  The  overall  sound  generation  processes 
have  been  classified  in  terms  of  the  sound  due  to  an  isolated 
turbulent  flame  and  that  due  to  the  interaction  of  a flame  with 
its  environment  in  a typical  combustor.  The  analysis  has  been 
focused  on,  (i)  the  far  field  noise  characteristics,  and  (ii) 
the  mechanism  of  sound  generation,  dispersion,  and  transmission 
in  the  vicinity  of  an  open  flame.  The  acoustic  intensity 
generated  by  a turbulent  premixed  flame  is  found  to  be  a 
function  of  the  relevant  aerothermochemical  parameters  and  the 
flame  structural  factor,  expressed  in  terms  of  six  double 
correlation  functions  characterizing  the  flame  structure. 
Explicit  expressions  for  the  sound  intensities  are  obtained 
based  on  a Wrinkled  flame  model  and  a Distributed  reaction 
model.  Noise  generated  by  liquid  droplets  are  classified  in 
terms  of  intrinsic  and  turbulent  driven  noise  components. 

The  intensity  of  the  intrinsic  noise  is  found  to  be  inversely 
proportional  to  the  fourth  power  of  the  mean  life  time  of  the 
droplet.  The  noise  amplification  by  acoustic  instability 
contributes  significantly  to  the  combustion  noise  in  high 
performance  ducted  spray  combustors. 

The  near  field  study  reveals  (i)  two  different  aerothermo- 
chemical roles  played  by  steady  and  non-steady  heat  release 
rate  with  regard  to  the  sound  generation,  and  (ii)  the  condi- 
tions for  the  resonant  oscillation  and  the  self-sustained 
oscillation  of  the  sound  wave. 
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NOMENCLATURE 
a = speed  of  sound 

A = nondimens ional  speed  of  sound  or  constant  defined 
by  Eq.  (Ill) 

An  = constant  in  the  solution  of  pressure  fluctuation 
B = constant  defined  by  Eq.  (Ill) 

1.1. 

= constant  in  the  frequency  factor  for  the  k1"  reaction 

Bn  = constant  in  the  solution  of  pressure  fluctuation 

c = constant  in  the  expression  of  the  ratio  of  turbulent 
burning  speed  to  that  of  laminar  burning  speed 

C = constant  defined  by  Eq.  (Ill) 

Cn  = constant  in  the  solution  of  pressure  fluctuation 
d = diameter  of  a droplet 

Dn  = constant  in  the  solution  of  pressure  fluctuation 

Ej,  = activation  energy  for  the  k*"*1  reaction 

f = non-dimensional  number  density  of  droplet 

fnj  = inhomogeneous  function  appears  in  Eq.  (20) 

P = wave  function  defined  by  Eq.  (33)  or  constant  defined 
by  Eq.  (112) 

3-  = function  defined  by  Eq.  (28) 

gnj  = inhomogeneous  function  appearing  in  Eq.  (19) 

G = constant  defined  by  Eq.  (112) 

h^  = specific  enthalpy  of  species  i 

hnj  = inhomogeneous  function  appearing  in  Eq.  (18) 

H = constant  defined  by  Eq.  (112) 

AH  = heat  of  reaction 

I = sound  intensity 

k = rate  constant 

K = function  defined  by  Eq.  (35)  or  kernel  defined  by 
Eq.  (93) 

Knj  = inhomogeneous  function  appearing  in  Eq.  (17) 
i = flame  thickness 

L = linear  differential  operator  defined  by  Eq.  (21)  or 

constant  defined  by  Eq.  (112),  or  latent  heat  of 
vaporization 

m = mass  of  a droplet 

m = mass  burning  rate  of  a droplet 
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velocity  component  in  transverse  direction 
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a,  = exponent  determining  the  temperature  dependence  of 
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1.  Introduction 


The  mechanisms  of  noise  generation  by  combustion  processes 
are  being  studied  in  depth,  since  a detailed  understanding  of 
these  processes  is  a key  to  providing  a description  of  the 
overall  noise  in  the  combustors  of  various  environmental, 
industrial,  and  aerospace  applications.  A qualitative  class- 
ification of  combustion  noise,  based  on  typical  frequency 
spectra,  reveals  random  noise,  registered  by  the  rough  burn- 
ing in  either  premixed  or  diffusion  flames,  superposed  on 
noise  at  selected  frequencies.  The  selected  frequencies  may 
be  caused  by  physico-chemical  interactions  occurring  between 
the  acoustic  field  and  the  system,  including  the  fuel  injec- 
tors, flame  holders,  and  the  flame  enclosure,  which  may  have 
been  excited  by  components  of  the  random  fluctuations  in  the 
flame.  While  the  selected  frequency  noise  can  be  eliminated, 
to  some  extent,  by  appropriate  designs  of  the  combustor 
system,  the  suppression  of  random  noise  appears  to  be  more 
complicated.  These  two  noise  fields  are  both  induced  by  ncn- 
steady  heat  release  in  the  burning  processes,  and  therefore, 
together  comprise  combustion  noise. 

Previous  studies  of  the  random  noise  generated  by  turbu- 
lent flames  have  been  made  by  several  investigators,  including 
Smith  and  Kilham  [1],  Hurle,  et  at  [2],  and  Giammar  and 
Putnam  [3] . Excellent  correlation  between  the  sound  pressure 
and  the  emission  intensity  of  C 2 and  CH  in  the  ethylene-air 
premixed  flames  [2]  demonstrates  the  intrinsic  nature  of  the 
turbulent  driven  random  noise  that  is  basically  different 
from  the  sound  generated  at  selected  frequencies.  Experimental 
studies  of  sound  generated  at  selected  frequency  bands, 
provoked  by  fuel  injectors,  or  air  entrainment  in  the  combus- 
tor, has  been  reported  by  Dance  and  Sutherland  [4],  Recent 
experimental  studies  with  a ducted  burner  conducted  by  the 
Princeton  group,  (Abdelhamid,  et  al,  [5]),  revealed  strong 
unsteady  combustion-cavity  coupling  resulting  in  intense  noise. 
Indeed,  the  noise  generation  by  the  interaction  of  non-steady 
combustion  with  its  enclosure  constitutes  a significant  pro- 
blem with  respect  to  core  engine  noise.  Basic  theory  on  the 
generation  and  amplification  of  noise  by  combustion  environ- 
mental interaction  is  well  summarized  by  the  books  published 
by  Gaydon  and  Wolf hard  [6] , and  Markstein  [7] . 

Examination  of  the  far  field  noise  patterns  has  established 
the  empirical  fact  the'_  the  source  of  combustion  noise  is  of 
monopole  type,  for  .joth  oremixed  and  diffusion  flames.  This 
observation  led  to  the  pnerw.aenological  theory  of  Bragg  [8] 
which  describes  the  flame  in  terms  of  monopole  sources.  The 
acoustic  power  and  thermo-acoustic  efficiencies  are  also 
obtained  from  phenomenological  reasoning,  with  resulting  good 
agreement  with  experimental  data. 
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Another  approach  to  the  problem  was  taken  by  Strahle  [9] 
who  adopted  Lighthill’s  formalism  for  a noise  field  predic- 
tion. The  first  decisively  important  factor  was  the  overall 
description  of  the  turbulent  flame  to  be  used  for  the  evalu- 
ation of  the  volume  integral  expressing  the  monopole  sound 
field.  Using  the  wrinkled  flame  model  developed  by  Karlovitz 
[10]  and  others,  Strahle  proposed  an  acoustic  power  scaling 
rule  similar  to  that  of  Bragg  but  differeing  in  two  ways, 
namely,  (i)  appropriate  pairing  of  the  exponents  appear  on  the 
gas  speed  and  the  flame  speed,  u4_33Sx)r+^<2,  and  (ii)  intro- 
duction of  the  scale  of  turbulence  aopropriately  powered, 

£2+r  , in  deference  to  the  phenomenological  wrinkled  flame 
model  which  has  large  scale  turbulence.  The  exponents,  r , 
and  q , appearing  in  the  scaling  rule  are  functions  of  the 
Reynolds  number  and  their  numerical  values  are  to  be  deter- 
mined from  the  experimental  data.  While  the  introduction  of 
these  empirical  factors  is  most  attractive  from  the  practical 
point  of  view,  the  underlying  physical  concept  of  the  correla- 
tion is  less  clear. 

The  scaling  rule  for  turbulent  flames,  based  on  the  model 
of  the  distributed  reaction  by  Summerfield  et  al  [11]  has 
also  been  used  by  Strahle.  There  is,  however,  reason  to 
believe  that  the  result  thus  obtained  has  limited  applicabil- 
ity because  of  simplifying  approximations  made  in  the  analysis. 

In  the  present  paper,  a general  formalistic  approach  to 
combustion  noise  is  developed  with  the  basic  physical  descrip- 
tions being  retained  at  a sophisticated  level,  and  the 
bifurcation  in  the  physical  as  well  as  mathematical  perceptor 
being  postponed  until  the  latest  possible  stage.  The  first 
part  of  the  paper  is  concerned  with  (1)  qualitative  and  quanti- 
tative descriptions  of  the  monopole  noise  sources  of  turbulent 
flames,  the  far  field  noise  characteristics  of  these  sources; 
and  (2)  qualitative  studies  of  the  generation  and  the  trans- 
mission of  sound  waves  in  the  near  field,  i.e.,  in  the  vicinity 
of  the  flame.  The  principal  result  of  the  present  investiga- 
tion of  flame  noise  is  that  an  explicit  far  field  noise 
expression  can  be  formulated  without  specifying  the  structure 
of  the  flame.  The  indispensable  feature  of  the  analytical 
step  is  the  derivation  of  the  wave  equation  for  the  reacting 
gas  in  the  most  general  terror-  The  wave  equation  obtained 
reveals  the  fact  that  the  heat  release  rate  contributes  to 
the  generation  of  sound  and  to  the  dispersion  of  the  sound 
wave.  Further,  on  investigating  the  effect  of  steady  and  non- 
steady heat  release,  it  is  evident  that  the  non-steady  and 
non-homogeneous  fluctuation  of  the  heat  release  in  the  turbu- 
lent flame  produces  six  dominant  sources  contributing  far 
field  acoustic  power.  These  six  sources  are  represented  by 
double  correlation  functions  obtained  from  the  fluctuating 
gas  velocity,  time  derivative,  and  spatial  derivatives  of  the 
fluctuating  heat  release  rate  together  with  the  weighting 
factors  comprised  of  steady  state  gas  velocity  and  the  spatial 
derivatives  of  heat  release  rate.  One  of  these  six  correla- 
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tion  functions  is  identical  to  the  acoustic  power  expression 
obtained  by  Strahle.  It  may  be  emphasized  than  the  sum  of 
these  correlation  functions  appears  expliciteiy  in  the 
scaling  rule,  and  is  termed  the  flame  structural  factor  in 
view  of  its  physical  significance.  This  factor  has  been 
used  in  the  scaling  rules  of  Bragg  and  Strahle  but  the  com- 
plete physical  implication  of  the  factor  in  reference  to  the 
flame  structure  has  not  been  recognized.  The  present  theory 
yields  an  intrinsic  structural  factor  which  may  be  evaluated 
by  an  experimental  survey  of  the  flame. 

While  the  general  characteristics  of  the  far  field  noise 
can  be  prescribed  by  the  correlation  functions  obtained  from 
the  formal  asymptotic  solution  of  the  wave  equation  for  the 
reacting  gas,  the  detailed  mechanisms  of  the  near  field  sound 
generation  and  dispersion  can  also  be  examined  within  the  flame 
and  its  vicinity  by  the  solution  of  the  wave  equation. 

The  near  field  problem  is  of  particular  interest  in 
reference  to  the  question  of  detailed  acoustic  phenomena  rela- 
tive to  imposing  turbulence,  the  aerothermochemical  structure 
of  the  flame,  effects  of  the  flame  holder,  or  upstream  condi- 
tion, and  the  condition  for  the  resonant  acoustic  oscillation. 

The  elucidation  of  these  phenomena  stems  from  the  consid- 
eration in  retrospect  to  the  potential  acoustic  roles  which 
non-steady,  non-homogeneous  heat  release  plays.  Indeed,  the 
self-contained  linearized  analysis  could  account  for  these 
intricate  phenomena  through  simplified  one-dimensional  steady 
flame  structure  superposed  by  three  dimensional  turbulent 
fluctuations.  The  near  field  analysis  has  been  conducted  on 
the  open  flame  of  known  flame  structure,  inclusive  of  steady 
and  turbulent  properties. 

While  some  research  has  been  aimed  at  obtaining  an  under- 
standing of  the  noise  generation  by  turbulent  flames,  little 
is  known  about  the  generation  of  noise  by  liquid  spray  combus- 
tion which  is  the  most  common  mode  of  combustion  in  commercial 
applications.  Liquid  spray  combustion  is  characterized  by 
various  complicated  non-steady  phenomena  occurring  concurrently 
in  a two  phase  system.  The  result  of  these  various  forms  of 
non-steadiness  is  the  fluctuation  in  heat  release  at  the 
frequency  of  the  dominant  mode  of  unsteadiness.  The  second 
part  of  this  paper  includes  results  of  an  investigation  of 
the  noise  generation  by  liquid  spray  combustion.  A detailed 
analysis  reveals  the  non-steadiness  in  the  combustion  of 
droplets  in  a time  scale  of  the  order  of  magnitude  of  the 
droplet  life  time,  generates  sound.  Furthermore,  the  turbu- 
lent fluctuation  alters  the  burning  rate  and  the  number 
density  of  the  droplets,  and  thus  generates  sound.  The 
acoustic  instability  in  the  liquid  spray  combustor  could 
yield  higher  rates  of  transfer  from  the  reacting  system  to  the 
acoustic  field  and  thereby  contributes  significant  sound  out 
put.  The  mechanism  of  sound  generation  originated  by 
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chemical  instability  is  of  particular  interest  in  view  of 
the  fact  that  the  fluctuating  reaction  is  the  basic  source 
of  sound.  The  existing  analysis  is  extended  to  consider 
the  effect  of  the  temperature  fluctuation  in  the  stability  of 
a homogeneous  reacting  gas  system.  The  analysis  yields  a 
basic  correlation  between  the  sound  intensity  and  the  chemi- 
cal parameters,  including  rate  constants  and  heats  of 
formation  of  reacting  species. 


2 Theoretical  Consideration  and  Mathematical  Formulation 
2.1  Theoretical  Consideration 

The  analytical  procedure  adopted  in  the  present  study  is 
aimed  at  two  basic  objectives  of  primary  interest  in  combus- 
tion generated  noise. 

The  first  objective  is  the  identification  and  quantita- 
tive description  of  noise  emitting  mechanisms  from  various 
combustion  processes,  including  homogeneous  and  non-homogene- 
ous  burning,  premixed  and  diffusion  flames,  open  or  enclosed 
flames.  The  salient  feature  common  to  these  various  noise 
emitting  mechanisms  and  noise  amplifying  combustion  processes 
is  that  the  sound  generation  is  attributed  to  the  non-steady 
intercoupling  between  various  aerothermochemical  modes  excited 
in  the  combustion  zone.  To  be  specific,  the  non-steady  exo- 
thermic and  endothermic  reactions  generate  random  noise,  and 
steady  combustion  in  the  presence  of  pressure  and  velocity 
fluctuations  result  in  noise  amplification.  It  is  evident 
that  a detailed  knowledge  of  the  intercoupling  between  vari- 
ous modes,  namely,  acoustic, entropy  and  vortex  modes  in  the 
combustion  zone  is  essential  in  the  identification  and  the 
quantitative  description  of  the  noise  sources.  For  this 
specific  reason  the  modal  study  developed  by  Chu  and  Kovasznay 
[12]  may  be  extended  to  the  burning  zone  characterized  by 
rapid  heat  release  and  gas  expansion. 

It  may  be  pointed  out  that  the  identification  and  the 
estimation  of  the  sound  sources  is  further  complicated  by 
interactions  of  the  flame  with  the  environment,  including  ducts, 
bluff  bodies,  fuel  injectors,  and  any  components  capable  of 
responding  acoustically  or  aerodynamically  with  the  combustion 
system.  The  subject  is  of  great  practical  Importance  in  the 
overall  assessment  of  engine  core  noise  on  which  active  re- 
search efforts  are  currently  underway.  The  analytical  pro- 
blems there  consist  of  sound  generation,  transmission  and 
reflection  in  an  environment  with  appropriate  boundary  condi- 
tions. General  mathematical  techniques  for  the  problem  of 
non -reacting  gases  were  developed  in  the  past,  however,  with 
the  exception  of  acoustic  instability  in  a rocket  combustion 
chamber,  little  effort  has  been  devoted  in  synthesizing  noise 
generation,  amplification,  and  transmission  in  ducted  combustors. 

The  second  objective  of  the  study  is  the  prediction  of 
the  far  field  intensity  of  the  sound  generated  by  combustion 
system.  The  analytical  schemes  for  the  open  flame  and  ducted 
flame  are  different  in  several  ways. 

The  evaluation  of  the  far  field  sound  intensity  of  open 
flames  is  complicated  by  the  non-isothermal  zone  and  the  non- 
isentropicity  prevailing  in  the  burning  zone.  In  view  of  these 
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two  physical  complexities/  Ligh thill's  theory  may  be  applied 
only  with  reservations.  An  alternative  approach  adopted  in 
the  present  analysis,  however,  is  the  method  of  multi-inter- 
vals, frequently  used  in  fluid  mechanical  problems.  Tvis 
method  offers  a self-consistent  procedure  of  predicting 
pressure  fluctuations  in  three  zones;  namely,  the  burning 
zone,  the  non-isothermal  layer  and  the  environment.  Solutions 
in  these  tl.ree  zones  are  matched  by  appropriate  techniques  for 
the  prevailing  physical  conditions.  The  mathematical  tech- 
nique for  the  prediction  of  the  far  field  noise  of  an  open 
flame  is  illustrated  in  section  3. 

Experimental  and  analytical  investigation  of  the  sound 
intensity  from  ducted  combustors  has  been  investigated  by  the 
Princeton  group.  The  experimental  correlation  of  the  internal 
pressure  fluctuation  and  the  far  field  sound  pressure  level 
agrees  well  with  the  theoretical  prediction  with  the  correc- 
tion accommodating  the  non-isothermal  environment. 

In  addition  to  the  objectives  stipulated  above, 

the  detailed  features  of  the  sources  of  noise  will  be  described 
in  an  effort  to  guide  the  analytical  study.  The  non-steady 
phenomena  responsible  for  noise  generation  are  largely  pro- 
voked by  various  mechanisms  of  fluid  mechanic  and  acoustic 
origin.  The  scale,  structure  and  intensities  are  different  for 
different  non-steady  phenomena.  Consequently,  various  sound 
fields  may  originate  from  the  same  combustion  processes.  For 
example,  droplet  combustion  is  inherently  unsteady  with  a 
characteristic  time  equal  to  the  order  of  the  life-time  of  the 
droplet.  The  combustion  of  a droplet  tr  quiescent  air  will 
therefore  generate  sound.  In  additior.  to  this  sound  field 
generated  by  a cloud  of  droplets , turbulence  with  an  eddy 
size  comparable  to  the  radius  of  the  flame  generates  noise  in 
much  the  same  manner  as  an  open  flame.  If  the  size  of  the 
.•>ddies  becomes  of  the  same  order  of  magnitude  as  the  inter- 
particle distance,  the  droplet  density  fluctuates.  This 
results  in  local  heat  release  fluctuations,  and  thereby 
sound  generation. 

The  last  observation  suggests  that  relevant  scaling  and 
non-dimensionalization  are  essential  in  identifying  and  screen- 
ing noise  sources.  This  will  be  illustrated  in  section  4.  All 
the  sound  fields  produced  may  be  adequately  superposed  for  a 
final  estimation  of  the  overall  sound,  provided  the  fluctuating 
pressure  field  is  substantially  lower  than  the  mean  pressure. 

If  the  fluctuation  becomes  of  a finite  strength,  the  super- 
position principle  ceases  to  apply  and  the  identification  of 
various  sound  fields  becomes  hopelessly  complicated.  Neverthe- 
less, the  perturbative  scheme  should  shed  light  on  the  physical 
mechanisms  of  the  generation  and  amplification  of  the  non- 
linear sound  field. 

It  has  been  recognized  that  many  reacting  systems  exhibit 
chemical  instability  by  which  the  rate  of  the  major  energy 
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releasing  reaction  fluctuates  at  characteristic  frequencies. 

If  the  resulting  fluctuation  in  the  heat  release  is  intense, 
the  sound  may  be  generated,  either  independently  or  coopera- 
tively with  other  non-steady  phenomena.  The  prediction  of 
the  sound  intensity  requires  a detailed  knowledge  of  the 
chemical  instability,  in  particular,  the  nature  of  the  heat 
release.  The  basic  understanding  of  the  noise  generation  by 
the  chemical  instability  is  also  deemed  as  a pertinent  object 
of  the  combustion  noise  theory. 

2.2  Mathematical  Formulation 

The  basic  equations  essential  for  the  study  of  combustion 
noise  are  the  conservation  laws  of  the  chemically  reacting  gas 
[13]  which  are  summarized  as  follows: 


2L  + 

pt 

y V-V  = o 

0) 

P d5  _ 
J Dt 

- Vf>  + V-X. 

(V 

5>  T = 

y 1 Dt 

$ + VKVT  - 

- jkcui/Wi)u>i 

(3) 

? Dyt  „ 
Dt 

•3 

> 

1 

.3 

3 

(4) 

where  S is  the  specific  entropy,  4>  is  the  rate  of  energy 
dissipation  per  unit  volume,  is  the  specific  chemical 

potential  of  i th  species,  wj[  is  the  rate  of  production 
given  by  the  phenomenological  chemical  kinetic  expression, 


V / w " 1i/  \o  -r**  -Et/er-fr (Xjf> 

e Jiff fj 

x vCtMi  - x vU'trbi 


where  w^  is  the  molecular  weight,  and  v"i  ^ are 

the  stoichiometric  coefficients  for  species  appearing  as  a 
reactant  and  product  respectively,  and  is  the  activation 

energy  of  the  kth  reaction.  It  is  assumed  that  the  gas 
mixture  obeys  the  ideal  gas  equation  of  state 

■P  = 

i = J 

While  the  conservation  laws  and  auxiliary  equations  are  self- 
contained  systems  for  the  present  study,  the  wave  equation 
which  describes  the  sound  field,  provides  mathematical  tools 
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for  a comprehensive  description  of  the  combustion  noise.  The 
derivation  of  the  wave  equation  is  elementary,  and  is  given 
in  Appendix  A.  The  resulting  equation  is  given  in  the 
following: 


dU.} 
* c>  Xj  53a 


f*T) 


Q « Zwthi  - 


•(8; 


The  remarkable  feature  of  the  above  equation  is  the  compact 
appearance  of  the  noise  generation  mechanisms,  including 
chemical  reaction,  heat  conduction,  diffusion,  and  viscous 
dissipation  which  are  the  processes  common  in  reacting  systems. 
(The  wave  equation  for  the  two  phase  system  is  given  in 
section  6 . ) The  second  term  appearing  on  the  right  hand  side 
of  Eq.  (7)  represents  the  source  of  the  combustion  noise. 
Special  attention  may  be  called  to  the  local  time  derivative 
and  the  spatial  derivative  of  an  effective  total  heat  release 
rate.  For  a flame  wherein  a rapid  rise  in  the  rate  of  the 
reaction  is  accompanied  by  expansion,  both  local  time  deriva- 
tives and  the  spatial  derivatives  of  the  rate  of  heat  release 
will  generate  a monopole  type  sound  field.  The  relative 
strength  of  these  two  terms  depends  to  a great  extent  on  the 
basic  flame  structure  and  the  turbulent  flow  field.  This  and 
other  aspects  will  be  discussed  quantitatively  in  the  section 
devoted  to  the  open  flame  noise. 

A qualitative  picture  may  be  drawn  with  regard  to  the 
difference  between  the  diffusion  flame  and  the  premixed  flame. 
Observation  of  the  sound  generation  source  reveals  that  the 
non-steady  diffusion  of  the  fuel  and  oxidizer  may  lead  to 
fluctuation  in  the  heat  liberation  and  the  location  of  the 
flame.  Furthermore,  the  non-steady  diffusion  may  contribute 
directly  in  the  quadrupole-like  noise  field  represented  by  the 
second  derivation  of  the  concentration.  The  order  of  magnitude 
of  the  noise  generated,  however,  may  be  small  compared  with 
the  combustion  noise. 
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3 Noise  Generation  by  Open  Turbulent  Flames 
3.1  Acoustic  Mode  in  the  Near  Field 

The  near  field  is  the  arena  dominated  by  strong  aero- 
thermochemical  coupling  between  a vortex  (turbulence)  mode, 
a sound  mode,  and  an  entropy  mode.  The  study  of  the  near 
field  will  be  focused  on  the  sound  mode  with  specific  refer- 
ence to  be  mde  to  the  turbulent  flame  zone.  This  suggests 
that  the  conservation  law  be  properly  scaled  by  the  tuibulent 
flame  thickness  if  . The  physical  variables  are  non- 
dimensionalized  by  appropriate  reference  quantities 

S;=  *i/*f , x = , a*=  a} /af  . M«=  \4/a,  , v£  = vc/v, 

where  , aj_  , are  the  velocity  and  the  speed  of  sound  at 
the  edge  of  tne  downstream  of  the  reacting  zone.  By  splitting 
the  velocity  V , into  components  consisting  of  basic  flow 
velocity  u and  perturbative_solencidal  component  ?<t>  , and 
rotational  field,  W , i.e.,  V = U + A<f>  + W , and  by  neglect- 
ing molecular  transport  processes  the  conservation  laws  and 
the  wave  equation  become 
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where  yi=-Mp;),A~-"x  p^uf  f,  ^His  the  heat  of  reaction  and  <^o 
is  the  rate  of  chemical  reaction  at  the  reference  state.  The 
non-dimensionalized  parameter,  A , represents  the  ratio  of 
the  heat  release  rate  to  the  thermal  convection,  i.e., 
A-^lV^/cpTi  Oi  where  is  the  burning  velocity.  Observation 
of  the  system  of  equations  Reveals  that  the  density  fluctua- 
tion at  low  Mach  number,  <<  1 , is  of  the  order  A , 
whereas  the  pressure  fluctuation  is  of  the  order  of  M^A  , i.e., 
the  compressibility  effect  associated  with  the  primary  density 
fluctuation.  This  suggests  that  the  perturbative  analysis  on 
the  basis  of  entropy  mode,  (A)  and  the  sound  mode  (M.^)  , 
should  render  a systematic  procedure  for  the  investigation  of 
the  near  field.  Accordingly,  the  flow  variables  are  expanded 
in  the  following  two-parameter  series: 

i(A.y'C)  (,4) 

J*0/J=0 
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The  velocity  field  is  expanded  likewise 


V = U + A ( V<p0il-t  Wo,,  ) + A2(V<l>0tZ+  Wo.i) 

+ tfZlA  (V<p2'l+  1*2.;)  + A*  (V<t>2.2  + ^2. 2)  + **  - Off; 


By  splitting  the  basic  flow  velocity  into  steady  and  non-steady 
parts,  i.e.,  U=u+u'  , the  total  derivative  becomes 


vtT  Pn  + Pt.  + Apr,.,*  A* vc0.i  * orvi"*'  %Vl. *)♦'* 


06) 


where  £ = Ww  + *..,)** 


An  important  physical  concept  pertaining  to  the  present  per- 
turbative scheme  is  that  the  small  Mach  number  perturbation 
scheme  is  valid  in  the  case  of  compact  noise  sources.  This  is 
evident  from  the  following  relation  Mo*  = V'/a,  - -ti /A  « 1 where 
X is  the  wave  length  of  the  sound.  The  noise  source  with 
Zf/X  <<  1 is  termed  as  an  acoustically  compact  source.  In 
general,  turbulent  open  flames  may  be  regarded  as  a compact 
noise  source. 

Substituting  all  the  variables  expressed  in  the  form  of 
Eq.  (14)  into  the  system  of  Eqs.  (9)  to  (13),  yeild  perturba- 
tive equations  for  the  sound  mode,  entropy  mode,  vorticity 
mode,  and  the  solenoidal  field. 


Sound  Mode  V*  4a  V/2n.j  = 

fin.j 

(17  ) 

Entropy  Mode  = 

h.  nj 

d? ) 

Vorticity  Mode  .E.Yffi  =. 

PT. 

3nj 

) 

Irrotational  Velocity  v*  j = 

(20 ) 

where  j = fnj /So  . The  inhomogeneous  terms  Knj  , hnj  , 

gnj  , and  fnj  are  given  in  Appendix  B.  * J 


The  wave  equation  governing  the  sound  mode  in  the  near 
field  is  degenerated  into  the  eliptic  system  in  this  region. 

The  fluctuation  which  drives  the  acoustic  field  acts  to  con- 
strain the  fluctuation  of  the  turbulent  field  nearly  dynami- 
cally incomprehensible  (m£  <<  1)  . The  inhomogeneous  term 
Knj  represents  various  mechanisms  including  true  noise  genera- 
tion, amplification,  and  refraction  pheno^na.  The  first  4 
terms  are  described  in  Table  1. 
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Some  relevant  phvsical  aspects  of  the  sound  mode  will  be 
described  categorically  as  follows.  Firstly,  the  primary 
sound  generation  appearing  in  the  M^A  and  M|»A2  approxima- 
tion is  attributed  to  the  turbulent  fluctuations  in  the  heat 
release  rate  and  the  gas  velocity  in  the  flame  structure. 
However,  the  fluctuations  may  also  be  induced  by  the  sound 
field,  and  thus  enhance  the  strength  of  the  noise  source. 

This  aspect  is  of  particular  interest  with  reference  to  acous- 
tic instability.  A term  Represented  by  Entropy-Basic  Flow- 
Sound  appearing  in  the  M“A*  approximation  is  the  potential 
source  responsible  for  the  chemico-acoustic  instability. 
Secondly,  the  sound  amplification  also  plays  a role  as  a re- 
sult ^f  entropy-acoustic  intercoupling.  The  maximum  amplifi- 
cation occurs  at  the  point  where  the  spatial  derivation  of  the 
heat  release  reaches  the  highest  value.  The  entropy  and 
vorticity  modes,  within  the  assumption  of  negligible  molecular 
transport  processes , are  convected  along  the  particle  path. 
These  modes  are  excited  primarily  by  heat  fluctuation  and  to 
a lesser  extent  by  sound  mode.  The  solenoidal  field,  repre- 
senting the  rate  of  volume  increase  per  unit  volume  of  the 
fluid,  is  also  primarily  affected  by  the  fluctuating  heat 
release. 

The  implication  of  the  observation  being  that  the  acous- 
tic mode  depends  on  other  three  modes  within  the  linearized 
approximation,  provided  the  pressure  fluctuation  remains 
small.  This  approximation  ceases  to  predict  a detailed  finite 
strength  sound  field,  particularly  when  combustion  instability 
dominates  the  burning  processes.  The  convergence  of  the  per- 
turbative scheme  also  has  to  be  established  whsn  it  is  applied 
in  the  non-linear  acoustic  regime. 

3.2  Intermediate  and  External  Zones 

3.2.1  Intermediate  Non-isothermal  Zone 

The  sound  waves  generated  or  amplified  in  the  flame  zone 
are  cransmitted  through  the  non-isothermal  layer  where  the 
waves  are  refracted  and  reflected.  In  the  absence  of  combus- 
tion, the  velocity  fluctuation  and  the  non-steady  heat  con- 
duction would  be  the  primary  sources  of  sound  generation  in 
the  non-isothermal  zone.  The  nature  of  the  sound  transmission 
depends  on  the  temperature  distribution  in  the  layer  and  the 
dimension  of  the  layer  relative  to  the  wave  length.  If  the 
characteristic  dimension  of  the  layer  exceeds  the  typical 
wave  length,  the  wave  emerges  from  the  high  temperature,  lower 
density  zone  and  propagates  into  the  zone  of  low  temperature, 
higher  density.  Physically,  the  non-isothermal  zone  may  act 
as  an  acoustic  cavity  with  a distributed  monopole  source  in 
the  inner  boundary  of  the  cavity,  i-o.,  turbulent  flame.  The 
waves  reflected  from  the  non-isothermal  layer  may  be  incident 
upon  the  flame  zone,  transmitted  through  it  and  reflected 
back  and  forth  inside  the  flame  zone,  (Figs.  1 and  2),  where 
the  temperature  is  substantially  lower  than  the  flame  zone. 
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1 schematic  diagram  of  sound  generation  by 
aflame : distributed  reaction  model. 
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Interaction  of  this  sound  wave  with  the  flame  may  provoke 
flame  oscillation  (Fig.  3a) . The  overall  effect  of  the  re- 
flected wave  may  not  be  ignored  for  a ducted  flame  with  pres- 
sure sensitive  exothermic  reactions. 

For  a small  size  burner,  the  overall  picture  is  substan- 
tially different  from  the  previous  case.  Within  the  limit  of 
the  thin  non-isothermal  layer,  (Fig.  3b) , the  overall  acous- 
tic effect  may  be  deemed  as  an  acoustically  compact  layer 
characterized  by  a rapid  decrease  in  the  gas  temperature. 

The  pressure  distribution  is  characterized  by  incompressi!  ie 
fluctuating  but  non-propagating  disturbances.  The  basic 
formulation  for  the  latter  problem  is  discussed  in  the  next 
subsection. 


Fig.  3a  Refraction  and 
reflection  of  sound  wave 
in  non-isothermal  acous- 
tic cavity,  X <<  Hq  . 


Fig.  3b  Refraction  and 
reflection  of  pressure  flue 
tuation  from  compact  non- 
isothermal  zone  X >>  5.0  . 


3.2.2  External  Zone 

A zone  extending  from  the  outer  edge  of  the  non-isothermal 
layer  to  the  geometric  and  acoustic  far  field  is  the  free 
space  for  sound  propagation.  The  primary  source,  if  it  exists, 
is  the  turbulent  velocity  fluctuation  due  to  the  jet  and  wake, 
which  is  neglected  in  this  study. 

While  the  solutions  and  the  matching  of  the  solutions 
for  these  three  fields  are  discussed  in  the  next  subsection, 
it  is  relevant  to  point  out  that  the  overall  problem  of  open 
flame  combustion  noise  may  be  regarded  as  consisting  of  a 
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thin  compact  noise  source  zone,  adjoined  by  a compact  non- 
isothermal  layer  or  by  a large  acoustic  cavity.  This  system 
is  finally  engulfed  by  an  extensive  wave  zone.  Matching  of 
the  solutions  from  a compact  zone  to  a wave  regime  calls  for 
the  singluar  perturbation  technique.  This  approach  is  self- 
consistent  and  should  provide  a detailed  aspect  of  the 
combustion  noise. 


3.3  Method  of  Matched  Solution  for  Compact  Flamelets 


The  equations  describing  the  fields  in  the  burning  zone, 
in  the  intermediate  non-isothermal  zone,  and  in  the  external 
zone,  will  be  solved  concurrently  and  matched  by  appropriate 
physical  conditions.  The  non-dimensionalization  of  the 
governing  equations  are  presented  in  Table  II,  and  the  re- 
sulting equations  within  the  first  order  approximation  for 
pressure  fluctuation  in  three  zones  are  given  as  follows: 


Inner  zone: 
L.JZ  = 


*5? 


-•/ 


htMtiM.  = gr&v 


•>! 


where  $<=Xi /Sf  , f.  = u;)4  , 


(2|) 


Intermediate 

A A 

L SI  ~ 


d/i  _ 


l i* 


6ru 


(22) 


where  ''l;  = Xc/Jle 


External  Zone: 


/■w 

L Si 


, c<~ 

4 X* 


36  6 


where  /K 


(23) 


The  basic  physical  conditions  for  matching  the  solutions  at 
zone  boundaries  are  the  continuity  of  pressure,  and  particle 
velocity. 

3.3.1  Solution  for  the  Inner  Zone 


The  general  solution  of  inhomogeneous  elliptic  Eq.  (21) 
may  be  solved  by  the  method  of  parametrix,  (see  for  example 
Courant  and  Hilbert  [14])  a method  developed  within  the 
formalism  of  the  integral  equation.  The  general  solution  is 
expressed  as 

Sl(  U,  x ) = 4 Jffi(SiiS!)JP<ti-sr)clV(SO  (J24-) 

where  is  the  homogeneous  solution  and 


V < f; : ?;'J  = -[(SrSiKSi-  SDJ’* 


is  the  parametrix. 
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Table  2 


Scaling  for  Turbulent  Flame  and  Spray  combustion  systems 


BURNING 

PROCESSES 


OPEN 
PREMIXED 
FLAME 
PISTRIBU 
7 GO  REAC 
VON  MODEL 


OPEN 
PR  E MIXED 
FLAME 
WRINKLED 
FLAME 

MO  PEL- 


OPEN 

SPRAY 

COMBUST 

LON 


OPEN 

SPRAY 

COMBUST 

LON 


DUCTED 

SPRAY 

combust 

LON 


TYPE  OP  FLAME 


SCALING  AND  FREQUENCIES 


INNER 

ZONE 


COMPACT 

SMALL 

SCALE 

TURB. 

L -4f 
uJ=“l 


SMALL 

SCALE 

TURB. 

L-U 

co~  Ut/Jf 

LARGE 

SCALE 

TURQ. 

L-'-Pt 

Ut  = Ui/f. 


INTERMEOt  EXTERNAL 
ATE  ZONE  ZONE 


WA  VE 
L s X 


WAVE 
L-  N 

U>=  a = J!l 


UL 

°y-.. 


• ••  . 
# • •«  * 


_ a _ -L. 

*°  “ X tUfc 


L-  j?t 

. Lh 


COMPACT 

L-Je 


WA  VE 
L -A 

W = £ s 

A Ja 


L-A 
u>~  ~ 


WAVE 
L = A 
u)~  a/j 


wave 

L = A 
to=  % 


TV  PE  OF 

SOU  NO 
FIELD 


MONOPOLE 


MONOPOLE 


uJ  =s  - 

X 

U 1 

compact 

WAVE 

L~Jq 

L-A 

uj=“± 

a , Ut 

■** 

< 

Mo  no POLE 


PL  POLE 


hiONOpOLE 


MONOpoLE 


MONOPOLE 


MONOpOLE 


MONOpoLE 


MONOPOLE 


The  unknown  function  ( 5 j,T)  has  to  be  chosen  so  that  the 
expression  for  J2.  given  by  Eg.  (24)  satisfies  Eg.  (21). 


In  view  of  the  specific  functional  form  assumed  by 
it  can  be  shown  that 

V2-/l  = 71/-  4-7T  o* 


£ / 


and 


L SL  - LJZ(a)- 


) (as*  b) 


Substituting  Eg.  (21)  into  Eg.  (25b),  yields  the  integral 
eguation  for  ) . 


YAMt>±  fOQ.  0$ ) 
4-7C  PZ  PZ  J 


U$) 


where  K(£^;|^')  is  given  by 


The  kernel  of  the  integral  eguation  has  a singularity,  and 
hence,  is  not  square  integrable.  However,  this  kernel  can  be 
reduced  to  a non-singular  one  by  a sufficient  number  of 
iterations.  Doubly  iterated  kernel  K, which  is  square 
integrable,  is  defined  as 


Kal J;:  ID  = & K ( t 

The  integral  equation  may  be  reformulated  in  terms  of 
kernel  as  follows: 


&( hrc)  =jff  K(zhL\  it)&( ii'c)  dv(tl) 


£3' 

- All^ 


The  equation  is  square  integrable,  and  thus  Fredholms'  theorem 
can  be  applied.  Since  the  exact  solution  of  Eq.  (27)  depends 
on  the  steady  state  gas  temperature  distribution,  no  attempt 
will  be  made  here  to  seek  the  solution  for  a specific  case. 
Instead,  a general  series  solution  may  be  written  formally  in 
the  following  form 


VA Mi]  / [do! 
4PC  M'U-c* 


provided  the  kernel  K(Ci?Ci)  converges  uniformly.  Kernel 
K(Ci;Cj)  is  the  sum  of  the  integral  operators 
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vXU; ; SO  = + /c%; si)  * Kulfc; si)  * ••* 


(aq) 


By  substituting  Eq.  (29)  into  Eq.  (24),  the  inner  solution  is 
given  by 


fU'Si.'L)  = J2!°'(f;,z)  + 


( 3o ) 


where  | /t-A'|  = C(fc-  $;)(  5;-  i:  3 ] /^2,  and  Ji  / ft  *.)  is 
the  second  term  which  appears  on  the  right  hand  side  of  Eq. 
(23) . It  is  interesting  to  note  that  the  third  term  appear- 
ing on  the  right  hand  side  of  Eq.  (3C)  is  of  the  order  of 


.-2. 


at  large  r.  In  fact,  in  the  geometric  far  field  within 


the  inner  zone,  Eq.  (30)  may  be  approximated  as 


J2(f-,x)  c:  71<s,(5o'c  > + 


y/i  rto* 


4 7C 


+ 0fO  (3I) 


The  primary  sound  field  is  a monopole  type  as  expected. 
3.3.2  Solutions  of  the  Non-Isothermal  Layer  and  Wave  Zones 


As  the  pressure  fluctuation  in  the  intermediate  zone 
depends  explicitely  on  the  temperature  distribution,  a quali- 
tative description  and  assumptions  with  regard  to  the  distri- 
bution of  temperature  will  be  given  prior  to  the  analytical 
investigation. 


The  temperature  distribution  in  the  vicinity  of  an  open 
flame  is  primarily  determined  by  convection,  conduction,  and 
turbulent  mixing.  While  convection  stretches  the  non-iso-- 
thermal  zone  in  an  axial  direction,  conduction  and  mixing 
transport  heat  in  a radial  direction.  Thus,  the  intermediate 
zone  may  be  approximated  as  nearly  spherical  in  shape  at  low 
convective  speeds.  This  approximation  is  in  accord  with  the 
unidirectional  characteristic  of  the  primary  sound  field. 
Departure  from  the  spherical  shape  of  this  zone  produces  a 
dipole  sound  characteristic  in  the  far  field. 


For  a spherical  intermediate  zone,  wherein  the  speed  of 
sound  is  a function  of  the  radial  coordinate,  the  solution 
of  Eq.  (22)  is  given  by 


J7.  CO.,'*-) 


K(t)  do. 


+ 


JT 


(32) 
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Where  n is  the  non-dimensional  radial  coordinate,  n = r/£g  / 
K(t),  CO)  M are  arbitrary  functions  which  will  be  deter- 
mined from  the  matching  conditions. 

The  solution  representing  the  outgoing  spherical  wave 
in  the  radiation  zone  is  given  by  the  expression 


XUJ.t) 


Ffr-  5) 

3 


where  £ is  the  non-dimensional  radial  coordinate. 


(33) 


3.3.3  Flame  Structural  Correlation  Functions  and  Far  Field 
Noise  Intensity 


Two  solutions,  Eqs.  (31)  and  (32),  are  matched  at  the 
interface  of  the  burning  zone  and  the  non-isothermal  zone. 
This  interface  is  located  sufficiently  far  downstream  of  the 
flame  where  the  chemical  reaction  ceases  to  exist,  but  where 
the  gas  temperature  is  approximately  equal  to  the  flame  tem- 
perature. It  will  be  further  assumed  that  the  interface  is 
located  in  the  acoustic  far  field  in  reference  to  the  reac- 
tion zone.  Let  the  location  of  the  interface  be  5 , in 
terms  of  inner  coordinate,  and  n , in  terms  of  the  inter- 
mediate coordinate,  i.e.. 

The  physical  boundary  conditions  are  the  continuity  of  the 
pressure  and  the  particle  velocity.  Application  of  these 
conditions  to  Eqs.  (31)  and  (32)  gives  the  following  expres- 
sions for  Kj[(t:)  and  3(0)  (^r) 


Jkfrta  M.J± 
Ki  ) Je 


osrj 


where  A^=  AQ  < Si)  . 


(36  ) 


Matching  Eqs.  (32)  and  (33)  is  carried  out  by  the  singular 
perturbation  technique  as  follows:  The  variables  rj  and 

are  related  so  that 


7 = 


(37) 


S 


Thus,  for  a fixedA  n , C 0 as  + 0 . ,Jhe  asymptotic 
expressions  for  ft(n,t)  with  n ■*  " / an^  n(?,r)  , for 
? ->-  0 , are  given  by 


n.  - 


7.  x)  . 


(38) 


and 

St  is.x)  ^ F<x)  + 4r1~>0 


(3q  J 
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The  continuity  conditions  give  the  follov/ing  relations 


(4e>) 


'XCtt  HA 

J2rt)  * 4^ 


fti) 


Substituting  Eq.  (41)  into  Eq.  (36)  gives 


-ft  i 

J1  C'Z) 


- * * Uhi  if'-  #*»] 

W%fffft£+ £)*««■>} 


Wet 

A2 


(42) 


The  far  field  solution,  in  the  wave  zone,  is  given  explicitely 
in  the  following  form 

*£■  * *$£(&* dVir!>  (43) 

where  6 ( i >-  T„ci) /Ttc  s *•  t,)  t Q '<■  Is.'Z)  is  the  temperature 

fluctuation.  This  expression  represents  the  monopole  sound 
field.  Comparing  the  far  field  expressions  in  the  radiation 
zone,  Eq.  (43),  with  the  burning  zone,  Eq.  (31)  , one  notes 
that  the  former  expression  is  similar  except  for  the  multi- 
plicative factor  of  the  ratio  of  the  square  of  the  speeds  of 
sound,  (A^/Aoo) 2 , or  Ti/Too  . The  appearance  of  the  tempera- 
ture factor  is  consistent  with  the  view  suggested  by  elemen- 
tary acoustics.  The  essential  similarity  preserved  in  far 
field  expressions  is  not  surprising  in  view  of  the  fact  that 
the  non-isothermal  layer  is  assumed  to  be  acoustically 
compact  so  that  the  refraction  and  reflection  of  the  pressure 
wave  in  the  thin  layer  does  not  alter  the  essential  emission 
characteristics.  This  suggests  a model  of  a composite  layer, 
consisting  of  a burning  zone  and  the  non-isothermal  layer, 
which  has  a source  strength  given  by  the  following  expression 

O-tfirAW*  f ] «**) 

As  the  thickness  of  the  non-isothermal  layer  increases,  the 
expression,  Eq.  (43),  deviates  from  the  actual  value  consider- 
ably in  view  of  the  refraction  and  reflection.  The  problem 
can  be  treated  by  replacing  Eq.  (22)  with  the  full  wave 
equation  which  can  be  obtained  by  appropriate  scaling  listed 
in  Table  II. 


The  far  field  noise  intensity  is  calculated  from  Eq. 
and  is  given  by  the  following  expression. 


(43)  , 
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- t* 

where  the  values  of  R^j 
double  correlation 


Rii  + %33  + *^a  + *^/3  + ^*1  ) 
are  obtained  from  the  following 


(AS) 


£ij  * Jf  Scj  (otjp ) dm)  dv(p)  lA 

VMV(f) 

a and  g denote  the  positions  of  two  correlation  elements. 
The  integrands  Sij  (a, 3)  are 


c*  _ aQw)  djk'tt)  c ,7 

~ n wr  » $«  = uLwu.(pj  jj. 


-i  - i.'. 


>33 


•s.  = , **  f'f",  s„  - * 

s<,=  i [ u-w- jf"’ 


The  intensity  of  the  far  field  sound  is  expressed  by  six 
sources  represented  by  double  correlation  functions  consisting 
of  the  fluctuating  velocity,  time  derivative,  and  spatial 
derivatives  of  the  fluctuating  heat  release  rate,  together 
with  appropriate  weighting  factors  comprised  of  steady  state 
gas  velocity  and  the  spatial  derivatives  of  heat  release  rate. 
In  view  of  the  physical  significance  attached  to  Rij  , the 
sum  of  these  functions  may  be  termed  as  the  flame  structural 
correlation  function.  The  factor  preceding  these  correlation 
functions  is  the  scaling  factor  for  the  noise  intensity.  While 
the  scaling  factor  depends  on  the  thermochemical  parameters 
and  the  upstream  condition,  the  structural  factor  contains 
flame  properties  which  depend  on  the  Reynolds  number,  (turbu- 
lence intensity),  and  other  aerothermo-chemical  parameters. 

The  evaluation  of  these  structural  factors,  experimentally  or 
analytically,  is  the  central  problem  of  the  open  flame  noise 
generation.  The  explicit  forms  of  structural  factor  based  on 
the  wrinkled  flame  model  and  the  distributed  reaction  model 
are  of  particular  interest.  This  is  described  in  the  follow- 
ing section. 
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4.  Structural  Correlation  Function  for  Wrinkled  Flame 
Model  and  Distributed  Reaction  Model 

The  contemporary  theories  of  turbulent  flames  are  well 
documented  in  many  books,  (see  for  example,  Williams  (13]  and 
Beers  (15]).  The  wrinkled  flame  model  advocated  by  Damkohler 
[16],  Schelkin  [17],  and  Xarlovitz  [10]  is  based  on  the  assump- 
tion that  the  flame  zone  is  a regime  in  which  the  fuel  is 
burned  as  in  laminar  flame.  Turbulence  is  considered  to  wrinkle 
and  perhaps  to  disconnect  the  flame;  and  augmentation  of  the 
burning  rate  results  from  the  enlargement  of  the  flame  area 
due  to  the  wrinkling  of  the  flame  surface.  The  distributed 
reaction  model  suggested  by  Summer fie Id,  et  al,  [11] , supports 
the  view  that  the  flame  structure  is  governed  by  the  small 
scale  turbulence  which  alters  the  effective  transport  coeffi- 
cient and  perhaps  the  effective  chemical  reaction  rate.  Struc- 
tural factors  for  these  two  flame  models  can  be  evciluated  in 
terms  of  the  turbulent  burning  speed  S>j«  . 

4.1  Wrinkled  Flame  Model 

The  generalized  view  of  the  wrinkled  flame  is  that  the 
flame  is  vnrinkled  by  large  scale  turbulence  and  the  flame 
structure  remains  essentially  laminar,  though  the  velocity 
and  temperature  may  fluctuate  within  the  flame,  as  in  Fig.  (4a) . 
An  element  of  non-steady  laminar  flame.  Fig.  (4b) , may  be 
considered  as  executing  linear  pulsation  and  rotation  by  the 
imposing  turbulent  fluctuation  about  a mean  flame  position. 

The  overall  wirnkled  flame  may  be  regarded  as  the  collection 
of  these  small  flamelets  with  their  translational  and  rota- 
tional motion  determined  by  the  local  turbulence.  The  total 
noise  intensity  is  calculated  by  the  sum  of  the  noise  pro- 
duced by  individual  filament. 

The  analysis  for  each  filament  is  carried  on  the  frame 
of  reference  O’  fixed  in  the  flame  structure;  (Fig.  4b) 
which  is  moving  relative  to  a fixed  point  0 . T simplify 
the  analysis  further,  the  rotational  motion  will  be  neglected. 

An  elementary  Galilean  transformation  yields  the  relation 
between  the  moving,  f , and  the  fixed,  x , coordinates. 

. __ 

X = I X(t  )ciz  + 5 = Xc(t)  + $ <48) 

~'o 

where  k(t)  is  the  velocity  of  the  point  C'  . All  the 
basic  scalar  conservation  laws  remain  invariant  under  Galilean 
transformation,  whereas  the  momentum  equation  becomes 

§£  = -fvPt  jrv-t  + m, 

The  v/ave  equation  is  found  to  be  invariant  under  this  trans- 
formation, and  thus,  the  results  obtained  in  previous  sections 
remain  applicable,  in  the  new  coordinate.  The  pressure  fluc- 
tuation is  the  far  field  is  given  by 
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Fig.  4a 

Model  of  Wrinkled 
Flame 


Fig.  4b 

Fluctuating 

Flamelet 


Fig.  4c 

Coordinate  System 


In  the  geometrically  far  field  |r|  >>  |x|  , |||  , Eq. 
(50)  is  expanded  in  the  following  form 


*n. 


_oq' 

~5x 


+ 


p'a 

PX 


)dmi) 


The  first  term  on  the  right  hand  side  represents  the  monopole 
sound  field  produced  by  the  small  scale  turbulence  in  the 
laminar  flame,  whereas  the  second  term  is  the  dipole  field 
generated  by  the  fluctuation  of  the  flamelet.  The  total 
noise  produced  by  the  wrinkled  flame  is 


-f--p  o MWo. 


-fo 


4-K 


i(£f{ip!f%£  *£)***'> 


c Vi 

If  the  flamelets  are  connected,  and  if  the  rate  of  the  reaction 
is  independent  of  time,  despite  the  velocity  fluctuation  (re- 
lative to  the  point  O'  fixed  in  the  flame),  the  monopole 
field  may  be  rewritten  as  follows 

- ygy&fcSfjt f dvt  w ‘ 5 1 > 

It  is  of  practical  interest  to  deduce  the  sound  intensity 
formula  in  terms  of  the  turbulent  burning  speed  Sip  . Note 
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that  Q 'v  uij?  = SL(de/d£n)  , where  c is  the  reaction  progress 
variable,  £n  is  the  coordinate  measured  in  the  direction 
normal  to  the  flame  surface.  Recalling  that  SjdA$  = S^dAy 
by  virtue  of  the  essential  hypothesis  of  the  wrinkled  flame 
model,  Eq.  (51)  reduces  to 


To  +*■ 


Thus,  the  sound  intensity  is  given  by 


where 


t Of- n(*H)jjsZvf  /sT  fpz  Cfc'3o. 

" /6  7Cxa£  A?-  \ SL  ) ' 

r=  J f'aJfl0;‘0Tn  (r3W 


and  dV  = dA^d^n  . Damkohler  and  Schelkin  have  derived  the 
ratio  of  the  effective  turbulent  burning  speed  to  that  of 
laminar  flame,  and  this  ratio  is  given  by  the  following  ex- 
pression , x 

For  high  intensity  turbulence,  the  following  expression  is 
used 

ST  = v' 

where  v'  is  the  turbulent  velocity  component  normal  to 
flame  surface.  The  ratio  of  turbulent  to  laminar  burning 
speed  is  experimentally  correlated  in  terms  of  Reynolds 
number  and  physico-chemical  parameters , see  Table  3 . 

By  substituting  Eq.  (54)  into  Eq.  (53a)  and  assuming 
that  the  velocity  fluctuation  v'  is  linearly  correlated  with 
the  Reynolds  number,  v'  = aRe  , the  expressions  for  the 
sound  intensity  becomes 

, I ( i ^^)px  (3=c{ft  C5€ft.) 

UKza&M  1 1 ISlJ 

For  high  intensity  turbulence,  the  sound  intensity  is 
given,  alternatively,  by  the  following  expression: 

£ («>>) 

J K>KzQ.fc  Ji*  1 Si- 


Other  expressions  of  sound  intensity  for  different  burning 
rate  formulas  are  listed  in  Table  3,  and  the  dependence  of 
the  intensity  is  illustrated  in  Fig.  5. 


Table  3 


Sound  Intensity  Based  on  Wrinkled  Flame  Model 


Investigator 

Burning  velocity 

a _ fir  /tn'oJ  Si*  . 

J b-\ ten ff.swr 

Damkohler 

Karlovitz 

St/s..-  1 Re  , low  vel. 

ST /SL-  ^ Rt=*Rt  , wl9h  vc). 

(i  + a(?e)z 

A*  n 1 

oc  Re 

Bollinger  & 
Williams 

sVs,=  o..*d°-zV~ 

d s burner  Tube  di&m  (Cm) 

,O.SL„  o.4t 
<5.324 -d  Re 

5y5u=0  + (3  Re)2,  » low  vcl* 

A + /3  Re 

Schelkin 

sVsL-  -fr-  R e=^f?4 , wl* 

Scurlock 

5t/Su=0  + W/lfi 

i + c*  y/& 

where  r2  is  given  by  Eg.  (53b). 


Fig.  5 Dependence  of  the  sound  intensity  on  Reynolds 
number  based  on  wrinkled  flame  model. 


4.2  Distributed  Reaction  Model 


For  small  scale  turbulence,  Summerfield,  et  al,  [11] , 
proposed  a "Similarity  hypothesis"  relating  to  the  Reynolds 
number  based  on  the  burning  speed,  flame  thickness,  and 

viscosity,  i.e. , „ . 

SrJj  = SuJl  s r(4>]  (S7) 

y-r  Vi.  ; 


where  is  the  effective  turbulent  diffusivity,  JtT  , 

are  the  thickness  of  the  turbulent  and  the  laminar  flame 
respectively,  and  f ( <f> ) is  the  function  of  the  fuel-air  ratio 
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for  a given  fuel.  Noting  Q=  the  effective  turbu- 

lent diffusivity  is  given  by  Prandtl  mixing  length  theory 

Vt  = Vl  ( 1 + S^V,v,/^)-^u(U-OiRcV')  {^B) 

where  e is  the  reaction  progress  variable  is  the  ef- 

fective mixing  length,  and  v'  is  the  turbulent  veolcity 
component  normal  to  the  flame  surface.  Substituting  Eqs.  (57a ,b) 
into  Ea.  (45)  gives  the  following  expression  for  the  noise 
intensity 

1 - hk'osjFI* Z W 

The  integrands  Sj_j  appearing  in  the  double  correlation 
function  Rj_j  are  given  bv 

S„  = 

S*  = a‘  { ff*w 

^33  = U‘-L  f UJ  “ 

S,2  = 

52,=  jaRe(Ui--|ui)i^';(v-|f;)3  Hi 


The  intensity  expression,  Ea.  (59)  , can  be  rearranged  in  the 
following  form 


J = A f az/?e  ( ft/  + fa  r Ru)  + Cfo  t- 

7Cl  o£a.1J}£  ^ 

t R33  ^ 


where  aReRi3  = Rl3  , a2Re2Rn  = Rn  , etc.  For  high  inten- 
sity turbulence,  the  source  terms  proportional  to  the  square 
of  the  Reynolds  number  dominate  over  two  other  terms. 


5.  Near  Field  Structure 

While  the  far  field  noise  intensity  is  described  by  the 
correlation  functions  obtained  by  a matched  asymptotic  solu- 
tion, the  detailed  mechanisms  of  the  near  field  sound  genera- 
tion and  dispersion  are  of  particular  interest.  The  pattern 
of  pressure  fluctuation  in  this  compact  zone,  unlike  that  of 
the  free  space,  is  coupled  to  the  turbulent  fluctuation.  The 
near  field  sound  is  comprised  of  primarily  pseudo-sound  field 
generated  by  the  fluctuating  heat  release.  This  pseudo-sound 
field  is  subsequently  joined  by  the  true  sound  field  in  the 
adjacent  zones.  The  near  field  analysis  of  the  open  flame  of 
known  structure,  inclusive  of  steady  and  turbulent  properties, 
is  described  in  this  section. 

To  facilitate  the  analysis,  the  domain  of  interest  has 
been  divided  into  three  zones;  non-reacting  upstream  region, 

I,  reaction  zone,  II,  and  hot  downstream  region,  III,  (Fig. 

6a ) • 


The  analysis  will  be  simplified  by  assuming  that  the 
temperature  and  the  velocities  are  uniform  in  the  three  zones, 
and  no  reaction  occurs  in  the  upstream  and  the  downstream 
zones  except  in  the  reaction  wave . The  turbulent  fluctuation 
is  assumed  to  be  one  dimensional  and  the  speed  of  sound  in 
the  upstream  zone  is  unity,  whereas  the  sound  speed  in  the 
reaction  and  downstream  zone  is  taken  to  be  CU  * a0  (6*  /©„)£ , 

The  equations  in  these  three  zones  are  given  in  the 
following: 


For  a 
- 00  < 


zone  extending  from  upstream  to  the  flame  front 
S < 0 

(61 J 


dVj., 


-i-  0 

MS  <*$»  ~ ° 


The  linearized  convected  wave  equation  in  the  flame  zone  is, 
0 < ? < 1 , 


Pt*  Mi  J* 


+ n<a(ij|^2+  siz 
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iM  N (if,  T.  J 


(6  2) 
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where  M2  = MoOf/Og)^  • 
of  the  flame  (1  < £ < °°) 


where  the  subscript,  2 , 
(II)  and  (III). 


The  wave  eqaution  in  the  downstream 
is 

(63  ) 

refers  to  the  quantities  in  zones 
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In  the  present  analysis,  the  upstream  and  the  downstream 
regions  are  assumed  to  be  of  semi-infinite  extent,  and  that 
the  sound  waves  generated  in  the  flame  propagate  upstream  and 
dov/nstream  without  reflections.  In  reality,  the  geometrical 
configuration  of  the  flame  is  such  that  the  upstream  region 
does  not  extend  to  infinity,  and  that  the  sound  originate 
from  various  parts  of  the  flame  may  incident  on  the  flame 
front,  as  is  illustrated  in  Pigs.  1 and  2. 


The  fluctuation  in  the  rate  of  heat  production  is  given 
by  the  following  expression 

The  solutions  in  the  upstream  and  the  downstream  regions  are 
given  by 

_ + -7—5  + SnJ 

-a.ci.x)  = S/4*,e 


C^h-c  - — 5 r £*) 

= Z 6>,  6 it  Hi  Vt. 


C6S-) 


The  constants  An  , Bn  , <5n  , and  en  are  to  be  determined 
by  matching  , ^3  with  &2  at  respective  boundaries. 


It  is  physically  evident  that  the  sound  waves  generated 
in  the  reaction  zone  will  propagate  in  the  upstream  and  dov/n- 
stream directions.  The  waves  will  be  partially  transmitted 
to  zones  1 and  3 , and  partially  reflected  at  the  boundaries. 
Hence,  the  general  solution  for  the  reaction  should  include 
the  particular  solution  as  well  as  tv/o  homogeneous  solutions. 
The  general  solution  is  given  by 
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(66) 


The  boundary  conditions,  at  5=0  and  £ = 1 , relate  to 
the  continuity  of  the  amplitude  and  the  gas  velocity.  Thus, 


at  5 = 1 , 

tv/o  conditions  give 

„ v S « 

A*  e = Cn  + 

and 

■ • a ( 

(61) 
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and  at  £ = 1 , continuity  conditions  give 


-1(7 fSV-eO 

= c„  e 
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+ Pn  e n * Fp„<0 
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. n ( ^Mz  \^Hi+MzV2  <„  a!!1  ta»  a«1)  djyo  r-7fl) 

t0"l~-r  Ml^Je  = cn>,ne  % pnAnc  + gf*  7 ' 


where  Fpn(l)  and  are  t*le  value  and  the  slope 

of  the  particular  integral  of  Eq.  (66),  evaluated  at  £ = 1 . 
From  Eqs.  (67),  (68),  (69),  and  (70),  Anei^n,  Bneien,  Cn  , 
and  Dn  are  calcualted  as  follows: 
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It  is  of  particular  interest  to  note  that  the  amplitude  of 
the  pressure  waves  are  all  proportional  to 
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and  that  the  contribution  of  the  second  tern  of  this  expres- 
sion diminishes  at  higher  frequencies.  To  examine  the  quali- 
tative aspect  of  this  noise  generation  term,  Fp(£)  , the 
explicit  expression  of  FD  will  be  calculated  when  the  source 
function  Nn(£)  is  given"  in  the  following  Fourier  series 

N d,"  )-  z. 

Substituting  the  above  expression  into  Ea. 

i SnTT* 


lar  integral  is  given  by 

Mi  K*  ( 6 


i(SnTCtt  ftt)  t74  j 

(66) , the  particu- 


^p*  = 


- I ) 


(1-  MtV?  )(L  SnK.-\*)(\,SnK  - 


Cri) 


It  may  be  noted  that,  within  the  present  approximation,  the 
inhomogeneous  solution  approaches  infinity  as  iSnU-’-Xn  d)  or 
Xn(2)  . This  resonant  oscillation,  however,  is  likely  to  be 
suppressed  by  the  dissipative  effects  which  are  not  considered 
in  the  present  analysis.  Some  of  the  typical  modes  of  outgoing 
pressure  waves  of  different  harmonics  (n  = 1 , 2)  induced  by 
the  first  few  modes  of  the  turbulence  Sn  = 1,  2,  3,  4,  and  5, 
are  sketched  in  Fig.  6b. 
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REACTION  ZONE 


Fig.  6b.  Generation  and  Propagation  of  Pressure  Waves 
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6.  Noise  Generation  by  Spray  Combustion 
6.1  Intrinsic  Noise  and  Turbulent  Driven  Noise 


The  combustion  of  droplets  in  the  combustor  has  been 
recognized  as  complicated  phenomena  invoked  by  heterogeneous 
aerothermochemical  systems,  with  mass  momentum  and  energy 
exchange  between  two  phases,  in  addition  to  the  exothermic  gas 
phase  reaction.  The  mechanisms  of  the  sound  generation  are 
also  highly  complicated  in  view  of  the  many  non-steady 
phenomena  concurrently  taking  place  in  the  combustion  processes. 
The  conservation  laws  and  the  derivation  of  the  wave  equation 
are  given  in  Appendix  C,  and  the  v/ave  equation  is  listed  in 
the  following 


{-£[(«; -W Ft  - CCJ/]  } 


> iHi  iHi  _ yJL  j-  jLla 
dxi  f aty 


ere ) 


where  m is  the  rate  of  the  droplet  vaporization  per  unit 
volume;  Fj_  is  the  force  per  unit  volume  exerted  on  the  gas 
phase;  Wj_  is  the  droplet  velocity  component;  Tg,  is  the  tem- 
perature of  the  liquid;  L is  the  latent  heat  of  the  vaporiza- 
tion; and  ip  is  the  rate  of  the  heat  exchange  between  two 
phases.  Two  mechanisms  of  noise  generation,  namely,  intrinsic 
and  turbulent  driven,  and  the  intensities  of  these  sound  fields, 
will  be  described  in  the  following  sections. 


6.2  Intrinsic  Noise  Generation 


Spray  combustion  is  inherently  non-steadv  in  the  time 
scale  equal  to  the  droplet's  life  time,  which  ranges  from  10“3 
to  1 sec  for  50  ^ 200  micron  sized  droplets.  The  major  physi- 
cal processes  responsible  for  the  noise  generation  includes, 
vaporization  accompanied  by  mass  and  heat  exchange,  followed 
by  diffusion  of  the  fuel  vapor,  and  the  non-steady  combustion 
in  the  gas  phase.  The  following  simplified  assumptions  will 
be  made  in  the  analysis  of  the  intrinsic  noise  generation. 
Firstly,  the  free  stream  flow  is  steady,  and  the  droplet  inter- 
actions are  neglected.  Secondly,  droplets  are  all  assumed  to 
be  spherical,  and  the  combustion  of  the  droplet  starts 
simultaneously  as  the  vaporization  of  the  droplet  starts.  In 
general,  time  lag  of  the  order  of  the  gas  phase  diffusion, 
ranging  10"®  to  10-2  sec,  exists  between  vaporization  and 
combustion.  Finally,  the  droplets  are  assumed  to  be  in  dyna- 
mic equilibrium  with  the  gas  phase,  and  the  transport  phenomena 
contributes  an  insignificant  effect  on  the  generation  of 
intrinsic  sound.  The  equations  are  non-dimensionalized  as 
follows : 


m!  £%  - M/(  f 1 
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CpTi 
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where  D/oz  - 

are  defined  as  Si  =v  xc/d 


and  the  non-dimensional  quantities 
o j "C  — ‘ A2"  ~ Qf/ = 
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w.  = Uirfdo/r.,;^),  Ai = (n0  yr>i  a h)/s,  crrt  . ? = 4/d0  > /»  n//v0 


The  reference  quantity,  do  , is  the  typical  droplet  diameter, 
is  the  life-time  of  the  droplet,  Nq  is  the  typical 
droplet  number  density  per  unit  volume,  mjj,  is  the  mass  of  a 
droplet,  and  a^  , pi  , are  the  speed  of  sound,  gas 
density,  and  temperature  of  the  hot  gas  in  the  burning  zone. 

The  effective  Mach  number  appearing  in  Eq.  (79)  is,  in  general, 
small,  and  thus,  the  first  term  which  appears  on  the  left 
hand  side  of  Eq.  (79)  is  negligible.  The  terms  which  appear 
on  the  right  hand  side  are  the  sound  source  terms,  including 
the  effects  of  mass  ejection,  vaporization,  and  combustion. 

In  estimating  the  far  field  intensity,  it  will  be  further 
assumed  that  the  characteristic  dimension  of  the  burning  zone 
and  the  non-isothermal  layer  are  all  small  compared  with  the 
wave  length,  and  that  the  temperature  distribution  in  the 
non-isothermal  zone  is  nearly  spherical.  Following  the 
analysis  similar  to  what  has  been  presented  in  Section  3 
the  far  field  intensitv  is  calculated  to  be 
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The  order  of  magnitude  of  this  intensity  of  the  monopole 
fiej.d  can  be  estimated  by  the  simplifying  assumption  that  the 
numerical  density  of  the  droplet  remain.';  constant,  and  that 
the  hot  gas  density  and  temperature  are  uniform.  The  change 
in  the  droplet  diameter  is  assumed  be  given  by  the  follow- 
ing experimental  formula  [13] 
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Substituting  Eq.  (81)  into  Eql  (80), 
finally  given  by 

_ (r-i)\MarniAH)(  B l n _ 

Max  6^ 7L x Ko  Q.oo/1 ' Qeo'tiifzJ 


the  sound  intensity  is 


where  V is  the  total  volume  of  the  combustion  zone,  and 
B is  the  characteristic  dimension  of  the  combustion  zone. 

It  is  interesting  to  note  that  the  intensity  is  inversely  pro- 
portional to  the  fourth  power  of  the  life  time  of  the  droplet. 
This  expression,  however,  is  the  maximum  intensity  that  would 
be  reached  when  all  the  droplets  burn  simultaneously  in  the 
period  of  the  droplet  life  time.  In  practice,  the  burning  of 
the  droplet  is  random  in  time  and  in  space,  thus,  the  intrin- 
sic noise  expression  should  be  weighed  by  the  appropriate 
statistical  factor  g(r).  The  final  expression,  Eq.  (82), 
will  be  accordingly  modified. 
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6.3  Turbulent  Driven  Noise  in  Spray  Combustion 


f ! 

i 


i 
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Turbulent  fluctuations  in  the  spray  combustion  zone 
produce  complicated  non-steady  combustion  processes.  It  is 
physically  evident  that  the  turbulence,  with  its  eddy  sizes 
smaller  than  typical  inter-droplet  distances,  will  affect  the 
number  density  of  the  droplets,  and  more  importantly,  the 
burning  rate.  The  generation  of  noise  by  turbulent  fluctua- 
tions can  be  treated  by  the  method  applied  in  the  preceding 
subsection,  except  that  scaling  and  non-dimensionalication 
need  to  be  modified  accordingly,  as  shown  in  Table  7. 


The  physical  dimension  is  the  size  of  til.,  turbulent 
eddy,  2-e  / and  the  reference  time  scale  is,  ?,e/ui  . The 
wave  equation,  Eq.  ( 78) , is  now  deduced  into  the  following 
form : 


■* 

is  the  reference 


burning  rate  of  the  droplet.  Rewriting  the  right  hand  side 
of  Eq.  (80)  in  terms  of  steady  and  fluctuating  quantities, 
the  far  field  expression  for  the  pressure  fluctuation  becomes 
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Tne  fluctuations  in  the  velocity,  the  density,  and  the  tempera- 
ture of  the  gas,  as  well  as  the  numerical  density  and  the 
burning  rate,  affect  the  strength  of  noise  sources..  Those 
quantities  are  either  to  be  determined  experimentally  or 
analytically.  If  the  characteristic  time  of  the  turbulence  is 
larger  than  the  mass  diffusion  and  the  heat  transfer  between 
the  gas  and  the  droplets,  the  quasi-steady  expression  of  the 
mass  burning  rate  may  be  used.  For  example,  the  burning  rate 
of  a droplet  with  convective  effect  is  given  by  the  following 
expression  [13] 


and  the  fluctuation  in  the  burning  rate  is  given  by 

. L > . 2KKdo  CpTi  6i /u  — (a  H )yoi°6  / 1 Vo  Pe  Li 
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At  higher  frequencies,  the  quasi-steady  approximation  is 
inapplicable,  and  therefore,  the  explicit  non-steady  burning 


5. 


— 


35 


rate  expression  should  be  used  in  estimating  the  high  fre- 
quency noise  characteristics. 

6.4  Acoustic  Amplification  in  Combustion  Zone 


In  all  the  problems  treated  in  previous  sections,  the 
sources  are  assumed  to  be  compact.  In  large  combustors, 
however,  the  size  of  the  combustion  zone  may  exceed  the  pre- 
vailing acoustic  wave  length.  The  heat  release  may  generate 
pressure  fluctuations,  as  well  as  amplify  the  sound  waves 
passing  through  the  combustion  zone.  The  detailed  picture 
of  the  generation  and  the  amplification  can  be  studied  for 
the  case  of  large  spray  combustors. 

Consider  a one-dimensional  spray  combustor  of  total 
length  Lc  , in  which  a combustion  zone  of  liquid  droplets 
extends  form  x = 0 to  L . The  chamber  is  closed  at  x = 0 , 
and  open  to  the  atmosphere  at  x = Lc  . Equation  (78)  is 
non-dimensxonalized  by  the  following  quantities  5;=x;/z_ 

-c  = ta,/u  , a*=  aVaf  , m,  = 

o*2  _ ) , v 6(q  da  dAjl)  » re  a<rj 

where  J2.  ^ijP-'PoV'Po  > e = ^ * 

The  effects  of  mass  addition  and  the  heat  of  vaporization  are 
assumed  to  be  negligible  compared  with  that  of  the  heat  re- 
lease due  to  the  reaction  which  is  assumed  by  the  following 
expression  • _ , ryL,0in- c fo,  , 

QC?,r;  =5  + zQ»n,t)  = 5^  + z-o 

The  steady  state  heat  release  is  constant,  and  the  fluctuating 
component  Sn(£)  is  assumed  to  be  known.  The  noi'-steady 
burning  rate  of  the  liquid  propellant  is  often  expressed  in 
terms  of  the  product  of  the  response  function  and  t.ie  pressure 
fluctuation,  as  for  example, 

QnH.t)  9 /And, t)S2, 

It  can  be  shown  that  the  analysis  of  the  generation  and  the 
amplification,  based  on  Eqs.  (86)  and  (87),  are  essentially 
similar  within  the  linearized  approximation.  In  order  to  re- 
tain the  generality  in  the  present  ~nalysis,  the  heat  release 
rate  given  by  expression  (86),  is  ar  apted.  The  solution  for 
the  pressure  fluctuation  is  assumed  i..  the  following  form: 

SLCS.'X)  = 2 <t>n('s)eLY)X’ 

Substituting  Eq.  (88)  into  Eq.  (85),  and  neglecting  the  terms 
containing  , the  following  equation  governing  <£n  , yields 

<t>n  t *n  <Pr>  = ^ Z a 

where  _ 

Ao  = fa1-  inreG[)/Af  , Q - i**  * = 
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The  distributed  burning  zone  extends  between  x = 0 and 
x = h . Equation  ( 89)  may  be  rewritten  in  the  matrix  form 

<p"  -f-  A2  + 0.N  jd  &o) 

where  <f>,  A,  N,  V , are  given  by 


( 1 0 ■ • -0  ■ 

1 

o 

o 
/ 

It 

.*/= 

0 z 0 

o A o 

. v 

[ 0 o • • n_ 

, o b 

f 

/ -qq  o 

••01 

r ° 

O • • • O 1 

K- 


0 2-iaQ  o 


V'  = 


A 0 
-<J2  A 


MIJ 


i 


0 


0 7i-  naQ  J 


l A J 


y 


The  general  solution  of  Eq.  (90)  consists  of  two  homogeneous 
solutions  <ptj  2.  / and  an  inhomogeneous  solution  V , which 

can  be  caJculated  from  the  follov;ing  integral  equation 

ctnd  the  kernel  K(f—  S' ) is  given  by 

j-A"1  l exp[-  6A(?-f'j]  («t3) 

where  A_/  is  the  inverse  of  matrix  A . 

The  matrix  integral  equation,  Eq.  (92),  is  Volterra's  type 
and  can  be  solved  if  the  inhomogeneous  term  is  absolutely 
integrable.  The  solution  of  Eq.  (92)  is  approximated  as 
follows  J 

i'(t)  = /(f)  + af  A/KCf-t'jV(S'Jfcyje/S' 
i l 6 t' 

*ajNK  Or-  Y)V<i‘)J  N K(s'-  f'V/ Ct")di"+  • • ft*) 

where  •/<»/=  (<ts') 
Note  that  the  convergence  of  the  series  solution,  Eq.  (94) , 
is  assumed  to  be  a given  parameter  a . Two  homogeneous 
solutions,  6i  and  <f>2  are  given  respectively  by  the 
following  exoressions 
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These  two  solutxons  satisfy  the  following  canonical  boundary 
conditions : 


Cp0)  (%-0)  = O 

<P‘Z>  = = i 


<p(ii  (5-0)  - 1 
<pU)t  (5-0)  ~ o 


The  general  solution  is  obtained  by  approximate  linear  com- 
binations of  the  inhomogeneous  solution,.  Eq.  (94),  and  the 
two  homogeneous  solutions,  Eqs.  (96)  and  (97). 

In  the  non-reacting  zone,  the  pressure  fluctuation  is 
governed  by  the  homogeneous  wave  equation,  i.e., 


tffi  _ A*-  jVl  _ 


_LC  - K ^ 5 ^ j 


The  solution  for  fi  is  given  by 

A'-'S.'C)  - ? )A  $ ;/3]e  l/v-c  C»o»] 

where  A and  B are  constant  vectors  to  be  determined  from 
the  boundary  conditions.  The  solutions  in  the  reacting  zone 
and  the  hot  non-reacting  zone  are  to  be  joined  at  the  inter- 
section of  the  two  zones.  The  boundary  conditions  are  given 
as  follows: 

j _ o u'  = o a doz) 

^ - n , XJ.'  — U!  <103> 


u =•  a 


t - O U - O ^ (10Z) 

5=  1 JX  = A > u'r,  U!  Cl05) 

S=  h 

Since  the  particle  velocity  vanishes  at  C = 0 , the  solution 
in  the  reacting  zone  is  given  by 

_/2.  ( * -c  ) = J c*>  A.  $ - ce?X$J  A A J J? x 

*<? 

. r / 

xi')V($')cx>A.yds'}C  + iVfj 

/■J  (105 

+ a/  NKcs-rjvmfw**'  + 

0 

+ * J^nm-vmvoj  NKn'-rjvuvfindr'T- ■ 

where  C is  an  arbitrary  vector,  to  be  determined  from  the 
boundary  conditions  at  £ = 1 . The  continuity  of  the  pres- 
sure at  £ = 1 is  given  by 

£c<W/\  - C&o\j  £40  AjOl/f  0/  Jfy  J c + -fa) 

. + a \nku-  t'J  VOOf  . 'WdJf'-t 

fi  r*' 

+ a2j  M/<(t-$oz($oJw(y-rjwM{(y)4SsC^M)A+faA;Q 


The  continuity  of  the  velocity  is 
A.  A ■+  A trw  A.  f A l/(  J A»*l  A.  J /Coo  A J ^ Jf  J“  C 

» •■  J JJ=I 


(107) 


r7l  m'J  (*NKCl-  r')V(t")fir')<*l=  (HCC0M.JA 

Jj=,  J- 

At  x = h , 3a'/,, j = O thus 

■+  (c-Oxi  JJ.  h.  ) O — o H°8) 


If  sUh  r*\TL  , <rr  (m-r  -£  ) 7Z.  , then  neither  sin  nor 

cos  h.  is  singluar.  Thus, 

a - - (cod  ah.  J \&nMh)  A 


(10^) 


From  Eqs. (106)  , (107)  , and  (10TV  the  vectors 

are  calculated  from  the  following  matrix  equations: 

FC  -t  H - [sinM-  (CosM.)(co6yUk)  -t  sin  Mb  J A 


\ (u o) 


J 


] 


G C t L = f M cos  A - A(sin  Cl)(cosMh)  Sin  Ah  } A 

Solving  for  A -:.nd  C,  there  yield 

A - LT'H  -tl(J~,F][WU,F-Gf'[i—WV~'H] 

B = -G°sMhTl(sinAh  )A 

C - [wu'f  -0]  [ L-  W u~‘h ] 

where 

r - Cos\-  CosAj  A ^/nXf'ccsAS'  V($')d$' 

o 

G •— -A sin\  -f  A s i n^J  A 1 sinAt'cos^X'Vd'jolt' 

0 

H = clJ^KCI- V)  V(V)-f(V)  cn'  -t 

L = “-iHj 

U = si  - (”00,5  f cos  '(sin  A h. ) 

W — MCos.U  - M.(Ai»M)(CC>sMh)(lSinMh)  ^ 

The  noise  generation  by  intense  acoustic  instability  is 
of  particular  importance  in  view  of  the  fact  that  the  ampli- 
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tude  of  the  wave  is  proportional  to  (Nrnt>  &H)/ cp  Tt  . Thi 

is  M|“2  times  greater  than  the  sound  generation  by  any  tur- 
bulent driven,  compact  noise  source.  For  a large  combustor 
with  a pressure  sensitive  reaction,  therefore,  acoustic  insta 
bility  should  be  adequately  controlled. 

The  mechanismr.  of  the  generation  of  the  higher  harmonics 
and  the  modulation  of  the  sound  wave,  as  illustrated  in  the 
analysis,  are  of  special  interest.  Since  the  velocity  and 
pressure  fluctuations  at  the  exit  plane  of  the  duct  affect 
the  nature  of  the  monopole-like  sound  field  for  a ducted 
burner,  it  is  of  practical  importance  to  design  an  optimal 
duct  minimizing  these  results.  Note  that  in  an  open  flame, 
the  acoustic  instability  is  not  likely  to  occur,  and  hence, 
the  noise  essentially  dominates  lower  frequency  regions. 


7.  Noise  Generation  by  Chemical  Instability 


It  has  been  demonstrated,  in  pre /ious  sections,  that 
the  monopole-like  sound  field  generated  by  a turbulent  flame 
is  primarily  attributed  to  the  small  scale  fluctuations  in 
the  rate  of  heat  release,  velocity,  and  temperature.  The 
estimation  of  fluctuation  in  the  rate  of  chemical  reaction 
is,  therefore,  essential  to  the  determination  of  the  source 
strength.  The  small  scale  in  the  turbulence  may  enter  the 
flame  from  the  upstream  region,  or  be  generated  within  rhe 
reaction  zone. 

The  stability  of  the  flame  structure  has  been  the  object 
of  intense  research  with  regard  to  turbulent  flame  theory. 
However,  the  question  of  chemical  stability  has  not  been 
investigated  to  a satisfactory  level  due  to  the  fact  that 
such  an  analysis  is  exceedingly  complicated,  and  that  the 
chemical  kinetics  involved  are  not  well  understood.  In 
the  analysis  of  chemical  stability,  (see  for  example  Bradley 
[18]  Glansdorff,  and  Prigogine  [19]),  the  question  has  been 
focused  on  the  chemical  stability  associated  with  the  response 
of  the  fluctuation  of  species  concentration,  whereas  the  effect 
of  the  temperature  is  ignored.  The  inclusion  of  the  temp- 
erature variation  is  essential  in  estimating  the  non-steady 
heat  release. 

Following  the  reaction  scheme  considered  by  Frank- 
Kamenetiky,  [20]  A + 3 

£ + N -*  0 + aN  ) 

A -<■  N B 

where  A represents  reactants,  B final  products,  and  M,  N 
intermediates.  For  a homogeneous  medium,  the  rate  of  reaction 
of  species  M and  N are  given  by 

P - -fc.cTjy*  Yh  - *tCT)yM  y* 

f = -*j<T)yMyN  - *sct)  y*yn  0«*b) 

Note  that  k^,  k , and  k3  are  rate  constants  which  depend  on 
activation  energy,  and  temperature  among  other  things.  The 
linearized  equation  describing  a small  departure  of  the 
concentrations  of  M and  N from  the  steady  state  including 
the  temperature  variation  are  given  in  the  following  non- 
dimensional  form, 


d Vjj 

ax 


- *3  Fa  V + a e ' 


dy  n __ 

dx  = 


k,  yAy^  + b&' 


(1150.) 


disk ) 


=s  yA  yM'  - yA  Y/l  t 0 ' U >£c.) 

Where  - 'Al/Ko  > ® - T'/To  > ^ = Kot/f>0 
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<*=  cpTo/t{  »/3  -(<Lhrt  bhi)/^  ? kq  is  the  numerical  va' ue  of, 
say,  at  reference  temperature.-  H is  the  overall  heat  of 
reaction,  h^_  and  h2  are  non-dimensional  heat  of  formation  of 
intermediates  M ana  N respectively.  The  frequencies  of  the 
system  are  given  by  the  following  characteristic  equations, 

u>3  - Xu)L  - r Ck.Ka  y/t  CbKjh,-  ak,hi.)y^/oi']  u>  -o  OU) 
where  y = (S/oC 

It  can  be  shown  that  when  a and  b are  set  to  zero,  the 
characteristic  equation  gives  two  pure  imaginary  roots 

u>  - ± i (K.Ki)*  yA 

Thus,  the  concentrations  of  M and  N fluctuate  at  the 
frequency  determined  by  the  reaction  rate  Kj_,  K3,  and  the 
steady  state  concentration  of  the  reactant. 

In  general,  three  characteristic  frequencies  are  found 
to  be  jl 

Ui, l3.  =1  ± i {k,k3  y/r  (bKih,  - atuHOZ/u  - ■£  ] * (m  ) 


60  3 » o 


Thus,  according  to  the  linearized  theory,  if  a and  b are 
sufficiently  small,  such  that  the  quantity  in  the  bracket 
remains  positive,  the  system  is  unstable  if  y>0,  or 


f 9K, 


Kj  + (K3  a/C* 

HTt  3T 


> o 


(/If  J 


The  Condition  stipulated  by  Eqn.(118)  suggests  that  the 
stability  of  the  reaction  depends,  to  a great  extent,  on  the 
derivative  of  the  rate  of  formation  of  the  intermediates 
M and  N.  This  stability  limit  of  this  reaction  scheme  is 
a function  of  gas  temperature,  among  other  things. 

The  instability,  or  oscillatory  behaviour  indicated  by 
the  reaction  schemes  typical  of  the  combustion  process, 
together  with  the  upstream  turbulence,  require  further 
investigation,  in  particular,  when  the  oscillation  becomes 
self-sustained  (19)  . 

The  source  of  noise  produced  by  the  chemical  instability 
induced  by  a reaction  scheme  suggested  by  £qn.(113)  can  be 
evaluated  from  Eqn.(7)  as  follows. 


O'-  off'  = <*-<)%{ *< Ya  % - vA  % + (#-  £ 


% Y»n<  + (-§j  £f*- -jf3J  9*  y„ nt)  i'j 
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The  order  of  the  magnitude  of  the  frequency  is  inversely 
proportional  to  the  reaction  line  which  is  po/k0  • The  first 
and  second  terms  appearing  in  the  curly  bracket  are  the  rate 
of  heat  release  associated  with  the  fluctuations  of  the 
intermediates  M and  N,  whereas  the  third  term  represents 
the  effect  of  the  temperature  fluctuation.  Essentially  the 
same  procedure  of  calculations  for  other  homogeneous  reactions 
may  be  used. 
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8.  Conclusion 


A theoretical  description  of  noise  generation  and 
amplification  by  combustion  systems  has  been  developed  on 
the  basis  of  theory  of  reacting  gas  dynamics  and  the  theory 
of  sound  generation.  The  principal  result  of  the  open  flame 
noise  generation  formulation  is  that  the  far  field  intensity 
of  the  sound  generated  by  acoustically  compact  open  flames 
can  be  predicted  by  the  flame  structural  factor  consisting 
of  six  double  correlation  functions.  The  explicit  forms 
of  the  structural  factors  are  obtained  for  two  flame  models; 
wrinkled  flame  model  and  distributed  reaction  model.  The 
experimental  determination  of  the  structural  factor  of  a 
flame  appears  to  be  of  particular  interest  in  the  prediction 
of  the  noise  intensity.  The  analysis  also  disclosed  that 
the  small  scale  turbulence  in  the  reaction  zone  is  the 
primary  source  of  the  monopole  sound  filled  whereas  the 
large  scale  turbulence  contributes  dipole  emission.  The 
burning  zone  is  dominated,  primarily  by  psudo  sound,  which 
acts  to  constrain  the  field  to  remain  dynamically  incom- 
pressible. The  rate  of  thermo-acoustic  energy  conversion  is 
determined  by  the  eddy  size  and  the  gas  velocity  among  other 
factors.  Consequently  a given  combustible  gas  burning  with 
intense  turbulence  of  small  eddie  size  is  likely  to  be  an 
intense  acoustic  emitter. 

Noise  generation  in  liquid  spray  combustion  process  is 
attributed  to  various  non-steady  phenomena  including  turbulent 
fluctuations,  intrinsic  non-steady  burning  characteristic 
of  droplets,  acoustic  instability  and  other  types  of 
instability.  The  excitation  of  the  specific  sound  mode  is 
largely  determined  by  the  interaction  of  the  burning  zone  and 
its  environment.  For  example,  liquid  spray  combusion  in  a 
ducted  burner  could  be  a potential  noise  emitter  when  the 
acoustic  instability,  or  resonance  phenomena  trigger  intense 
fluctuations  in  the  heat  release.  The  investigation  of  such 
non-steady  burning  processes  in  a duct  is  essential  to  the 
overall  assesment  of  core  engine  noise. 

Noise  may  be  generated  independently  by  chemical  in- 
stability or  cooperatively  by  various  instability  phenomena. 
The  effect  of  the  chemical  instability  on  noise  generation 
is,  to  some  extent,  obscured  by  the  uncertainty  in  the 
chemical  kinetics  of  the  major  energy  yielding  reactions. 
Further  research  is  necessary  to  assess  the  over  all  noise 
contribution  by  chemical  instability. 
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Appendix  A 

Convected  Wave  Equation  of  a Reacting  Gas 


The  convected  wave  equation  of  the  reacting  gas  is  derived 
from  the  conservation  laws  as  follows.  The  overall  continuity 
equation  for  a reacting  gas  is 


where 


SL  + f iu>  - 0 

Ot  * r 6 AC  ° 
PX 


(A->  ) 


is  the  pressure,  S is  the 


_d 

and  p = p(p,S,Yj[)  where 

entropy  and  Y^  is  the  concentration  of  species 
momentum  equation  for  the  reacting  gas  is  given  by 

L J>P  . X 

Ob  f d XC  S d Xj 

v;here  Tij  is  the  viscous  stress  tensor. 


ft  ^ ft 


The 


{A  -Z 


To  obtain  the  convected  wave  equation  for  the  reacting 
gases,  the  divergence  of  the  momentum  equation  and  substantial 
derivative  of  the  continuitv  eauation  is  taken  and  combined. 

THatiS'  i fflL.-l._i.fj.it  . 

and  expanding  the  density  derivative  in  terms  of  its  func- 
tional variables,  the  continuity  equation  becomes  after  some 
reduction  , . ^ Q ^ 

- i &*(#)**& 

By  subtracting  equation  (A-3)  from  (A-4)  and  after  applying 
some  rearrangement,  the  convected  wave  equation  for  a react- 
ing gas  is  obtained  in  the  form 


For  a reacting  gas,  Ec. 
ing  form 


(A-5)  is  reduced  to  the  follow- 


where 


Q = Z w;  hi  - 2 j;,  hi,  - Cl/;  - Vj  4\» 


V]_  is  the . diffusion  velocity  of  species  i,  Dx  • is  the 
thermal  diffusion  coefficient  for  species  i,  and  is  the 
viscous  dissipation.  * v 
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Appendix  B 


Perturbation  Equations  for  Acoustic  Mode, 

Entropy  Mode,  and  Vortex  Mode  in  the  Burning  Zone 

The  perturbation  equations  for  the  burning  zone  are 
obtained  by  substituting  the  assumed  tv/o-parameter  solutions 
of  the  flow  variables,  (14),  into  Eqs.  (9),  (10),  (11),  (12), 
and  (13).  By  collecting  the  terms  associated  with  the  same 
power  of  Moo  and  A , they  yield  systematic  perturbation 
equations . 

Acoustic  Mode: 

V-/4oVJ2n.j-  = Kr>.j 

J<z,i  = -JT 

l ^ r & • 


Entropy  Mode: 


= h-i 

flo.1  = - « - 

ho.  = - a.„-  &o*  A..;  - tyiii 

►>*•■  = y(  £,5'v.Ki  - ) 

h*’1-  r (^+a^tt«+  {®  M-Aotati.')} 


Rotational  Mode: 
(D  * n'.  n 
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Dilation: 


,£>jl  “ */o., 


V/22.i  ~ (£*uV)\74>0. 


0 fo.O 

9a.i  =.  --y-  vyi4.(-(^tu'7)74»2., 

0 fOiO 

92.2  = - TT-  ( ™4-2  + 2 ■««  - °Z.i  VA,.|  - <n.,  VJlz.,) 

a Jo.o  ' 


Field 
2 


V 4>».j  = 

1>,1=  2'  , = 

f0.2  = Qo.l  , f2.2=-J-(^l2,'t  T 82.1 
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Appendix  C 

Wave  Equation  for  Liquid  Fuel  Spray 


For  a two-phase  system,  the  conservation  laws  given  by 
Eqs.  (1)  - (4),  need  be  supplemented  by  the  conservation 
laws  for  a condensed  phase  and  a set  of  equations  describing 
the  rate  of  exchange  of  mass,  momentum,  and  energy  between 
two  phases,  see  for  example  Sirignano  [21] . A straight- 
forward algebraic  manipulation  shows  that  three  additional 
terms  enter  into  the  right  hand  side  of  Eq.  (A-5) . These 
three  terms  represent  the  effects  of  mass  addition, 
distributed  body  force,  -f-fccOVf)  . and  the  entropy 
generation  m'lr.p.  due  to  non-equilibrium  phenomena,  and 
exchange  of,  thermal  energy  between  two  phases.  The  rate  of 
entropy  •pCp£)t*p.  is  given  by  the  following  expression 


s P, f/ir-wl  DiU£-  a 

Pt JT.p.  "'JcTT  x dTJ* 


ovx  L 
fTi 


Cc-i  J 


Addition  of  (C-l)  in  the  right  hand  side  of  (A-5),  which 
contains  the  entropy  generation  due  to  gas  phase  reaction, 
gives  the  wave  equation  for  liquid  fuel  spray  combustion. 
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